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ABSTRACT 


A  new  concept  for  the  design  of  aeroservoelastic 
systems  is  introduced:  flutter  taming  by  nonlinear  control* 
i.  e.  *  use  of  nonlinear  terms  in  the  equation  to  ensure  that 
the  behavior  of  the  system  beyond  the  flutter  speed  is  of 
benign  rather  than  destructive  nature.  This  is  accomplished 
by  using  a  very  simple  nonlinear  control  law.  It  is  shown 
(using  a  singular  perturbation  analysis  about  the  stability 
boundary)  that  flutter  taming  is  always  possible  for  an 
aeroservoelastic  system  tliat  can  be  represented  by  a  system 
of  nonlinear  differential  equations  with  analytical 
nonlinearities.  It  is  important  to  emphasise  that  the 
control  system  for  flutter  taming  is  fully  nonlinear*  and 
therefore  it  does  not  affect  the  linear  behavior  (in 
particular  the  stability  characteristics)  of  the  system. 
Hence*  flutter  taming  can  be  used  in  conjunction  with 
flutter  suppression  by  active  control  (i.e. *  use  of  linear 
active  control  to  increase  the  flutter  speed).  Applications 
of  the  theory  to  the  case  of  an  airfoil  in  supersonic  flow 
are  presented.  In  addition  to  an  active  control  modification 
(use  of  control  surface  with-  nonlinear  feedback)*  passive 
modifications  (e. g. .  a  'nonlinear  damper)  are  also 
Investigated. 
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SECTION  X 
INTRODUCTION 


The  linear  description  of  flutter  and  the  use  of  linear 
active  control  to  increase  the  flutter  speed  (flutter 
suppression)  are  assumed  to  be. known  and  therefore  neither  a 
description  of  the  phenomenon  nor  a  review  of  the  extensive 
literature  in  the  field  are  provided  here  (see*  for 
instance*  Refs.  1*2  and  3).  Rather*  the  phenomenon  of 
interest  here  is  the  behavior  of  an  aeroservoelastic  system 
beyond  the  flutter  speed*  i. e. *  in  the  range  for  which  the 
linear  analysis  predicts  unstable  behavior.  The  nature  of 
nonlinear  flutter  was  considerably  clarified  by  the  research 
on  nonlinear  panel  flutter  which  includes  the  pioneering 
works  of  Librescu  (Ref.  4)  and  Kobayashi  (Ref.  9)  as  well  as 
the  extensive  work  of  Dowell  (Ref.  6  -  10).  The  results  of 
this  research  activity  indicate  that  the  post-flutter 
behavior  may  be  characterized  as  follows:  the  amplitude  of 
the  oscillation  of  the  panel  beyond  the  flutter  speed  grows 
in  time  (in  agreement  with  the  instability  predicted  by  the 
linear  analysis)  but*  in  contrast  to  the  linear  analysis 
prediction*  the  amplitude  growth  often  tends  to  a  finite 
valuej  the  limiting  value  solution  is  called  limit  cycle. 
The  above  phenomenon  is  called  benign  flutter.  Zf  the 
amplitude  grows  beyond  any  limit*  the  phenomenon  is  called 
destructive  flutter.  The  mechanism  of  nonlinear  panel 
flutter  was  clarified  from  a  mathematical  point  of  view  by 
this  author  in  Ref.  11.  Because  of  the  considerable  bearing 
on  the  results  presented  here*  the  above  work  and  its 
extensions  are  extensively  reviewed  here.  In  Ref.  11*  an 
analytic  expression  of  the  transient  response  was  obtained 
in  the  form  of  an  asymptotic  solution  using  a  singular 
perturbation  method  (known  as  multiple  time-scaling  method) 
which  had  been  developed*  apparently  independently*  by 
Sandri  (Ref.  12)«  Cole  and  Kevorkian  (Ref.  13)*  and  Nayfeh 
(Ref.  14). 

The  results  indicate  that  in  the  neighborhood  of  the 
instability  boundary*  an  approximate  solution  (with  the 
order  of  magnitude  of  the  error  equal  to  the  fourth  power 
of  the  limit  cycle  amplitude)  is  obtained  for  the  transient 
response  of  the  paneli  the  nature  of  this  solution  is 
characterized  by  one  single  real  number  (indicated  with  ^ 
in  Ref.  11).  If  such  number  is  positive*  then  the  nonlinear 
terms  have  a  stabilizing  effect  and  they  counteract  the 
effect  of  the  linear  terms  (which  are  destabilizing  above 
the  flutter  speed)  to  yield  benign  flutter  (stable  limit 
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cycle).  On  the  other  hand*  if  such  a  number  is  negative* 
then  destructive  flutter  occurs:  it  should  be  emphasized 
that  in  this  case  even  for  speeds  of  flight  below  the 
flutter  speed* where  the  linear  analysis  predicts  stable 
behavior*  the  amplitude  may  grow  beyond  any  limit  provided 
that  the  initial  amplitude  is  "sufficiently"  high: from  a 
mathematical  point  of  view  this  behavior  is  characterized  by 
the  presence  of  an  unstable  limit  cycle.  The  mathematical 
description  of  nonlinear  flutter  is  presented  in  detail 
later  in  this  paper. 


It  should  be  emphasized  that  the  results  of  Ref.  11  are 
not  limited  to  panel  flutter:  as  shown  in  Ref.  15.  the 
results  apply  to  any  dynamic  system  which  can  be  represented 
by  a  system  of  nonlinear  differential  equations  with 
nonlinear  terms  of  analytic  nature.  In  addition.  the 
analysis  can  be  extended  to  include  the  effect  of  fifth 
order  nonlinearities  (Ref.  16).  Finally.  the  same  results 
are  obtained  by  using  a  completely  independent  technique, 
the  Lie  transform  method  (Ref.  17). 


Two  slightly  different  formulation  of  flutter  analysis 
are  used  in  Refs.  11  and  15.  respectively.  The  first  one. 
presented  in  Ref.  11  deals  with  dynamic  systems  which  can  be 
represented  by  a  system  of  second  order  differential 
equations.  This  formulation  (presented  in  terms  used  for 
structural  dynamics  problems)  is  convenient  to  study  for 
instance  flutter  taming  of  an  aircraft  that  is  not  subject 
to  linear  active  control  (in  particular*  without  flutter 
suppression):  such  a  system  will  be  referred  to  as  an 
aeroelastic  system.  in  contrast  to  a  aeroservoelastic 
system,  which  includes  active  control.  Since  the  equations 
governing  active  control  are  usually  given  as  a  set  of  first 
order  differential  equations  (state-variable  approach),  in 
order  to  study  flutter  and/or  flutter  taming  of  an  aircraft 
subject  to  active  control.  it  is  convenient  to  recast  the 
mathematical  model  as  a  set  of  first  order  differential 
equations.  This  system  is  analyzed  in  Ref.  15  (where  the 
equivalence  of  the  two  formulations  is  also  discussed).  This 
second  approach  is  used  here  because  of  its  generality  and 
its  formal  simplicity. 


The  objective  of  this  work  is  to  show  (via  a  singular- 
perturbation  analysis)  that  nonlinear  feedback  control  may 
be  used  to  "tame"  the  phenomenon  of  flutter,  that  is  to 
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ensure  that  only  benign  flutter  (stable  limit  cycle)  occur 
and  that  the  limit-cycle  amplitude  be  below  any  prescribe 
value.  The  general  theory  of  Ref.  19  is  briefly  reviewed  in 
Section  2  and  used  to  discuss  flutter  taming  in  Section  3. 
The  methodology  is  then  applied  in  Section  4  to  the 
particular  case  of  an  airfoil  in  supersonic  flow.  Numerical 
results*  presented  in  Section  9*  indicate  that  acceptable 
values  for  the  limit-cycle  amplitude  may  be  obtained. 
Concluding  remarks  are  presented  in  Section  6. 


SECTION  II 


NONLINEAR  FLUTTER 


5L_1  Introduction 


Consider  an  aeroservoelastic  system  which  may  ba 
described  by  a  aat  of  aquations  of  tha  typa 

A*  /  dk  =  A  (A)x  +  f  (x,  \)  ( 1 ) 

whara  x  is  tha  vactor  of  tha  stata  variablas  tfacribing  tha 
systems  (a. g. »  generalized  Lagrangian  coordinates  and 
velocities  as  well  as  control  variablas).  In  addition  A(A) 
is  a  square  matrixi  A  is  a  positive  parameter  (a. g. ,  tha 

spaed  of  flight  or  tha  dynamic  pressure)  and  f  is  tha  vactor 

comprising  all  tha  nonlinear  terms  which  is  assumed  to  ba  an 
analytic  function  of  x  so  that 

f  s  j  £  VA)  *<■  xf  +£  v  *1  *'*■"}  <a’ 

The  assumption  under  which  Eq.  1  provides  a  satisfactory 
representation  of  a  typical  aeroelastic  system  are  examined 
in  Appendix  A.  For  very  small  amplitudes  the  system  is 
approximated  by  the  linearized  equation 

d  X  I dk.  -  }  X  (3) 

It  is  assumed  that,  for  A»  0,  the  system  is  stable,  in 
the  sense  that  all  the  eigenvalues,  s,  of  the  matrix  A(0) 
have  a  negative  real  part.  It  is  also  assumed  that  as 

A  increases,  at  A«  A0  one  (and  only  one)  pair  of  complex 

conjugate  eigenvalues  of  A<\)  cross  the  imaginary  axis.  Let 

(4) 

be  such  an  eigenvalue:  <*Jt  is  called  the  flutter  frequency 
whereas  A0  is  called  the  flutter  boundary  or  stability 
boundary.  Summarizing,  we  may  say  that  for  0  £A<A,  all  the 
eigenvalues  have  a  negative  real  part.  For  A#<  A <  A'  (where 
A'  is  the  value  of  A  at  which  another  crossing  of  the 
imaginary  axis  occurs),  one  pair  of  complex  conjugate 
eigenvalues  has  a  positive  real  part  and  the  others  have  a 
negative  real  part.  For  A  ■  A0  one  pair  of  complex  cojugate 
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parts.  Recall  that  (assuming  for  simplicity  that  tha 
eigenvalues  are  distinct)  the  solution  of  the  linear  system 
is  of  the  type 

r-i  5;  h 

X  -  2-  C.  £  U.  (9) 

—  ^  X  -  ' 

inhere  c^  are  constants  (which  depend  upon  the  initial 
conditions)*  whereas  s^  and  are  the  eigenvalues  and  the 

corresponding  eigenvectors  of  the  matrix  A(X).  that  is.  by 
definition,  are  such  that 


[ 5.  I  -  /V\)J  U-  r  0 


(6) 


with  ly^  I  #  0.  Then  we  may  say  that  for  0<A<Ae  the  solution 
of  the'linear  system  is  stable.  for  ~  A#<  X  <  A/  it  is 
unstable,  whereas  for  A*  A,  the  solution  is 'at  the  stability 
boundary:  as  time  goes  to  infinity,  the  contribution  of  the 
negative-real-part  eigenvalues  goes  to  zero  and  therefore 
the  steady  state  solution  is 


l  tvj  t  «  ^  t  V 

x  «  fr  e  •  +  c*  e  u 


(7) 


where  #  indicates  the  complex  conjugate  value.  c  is  a 
constant  which  depends  upon  the  initial 
the  imaginary  eigenvalue  and  u  is 

eigenvector  so  that  ” 


conditions.  iu>m  is 
the  corresponding 


[•iH  I  -  A.J  w.  --  2 


with 


(B) 


A.  --  A  fx.) 


(9) 


In  essence.  the  multiple-time-scale  analysis  used  in 
Refs.  11  and  19  consists  in  finding  the  solution  for  values 
of  A  in  the  neighborhood  of  A,,  using  a  singular  perturbation 
of  the  steady-state  solution  at  A*A,  (given  by  Eq.  7).  This 
analysis  is  presented  in  simplified  form  in  the  following 
subsections. 


Two  cases  are  not  included  here.  The  first  is  the  case 
in  which  at  A«A#  a  real  eigenvalue  (instead  of  a  pair  of 
complex  conjugate  eigenvalues)  crosses  the  imaginary  axis. 
In  aeroelasticity.  this  instability  is  known  as  divergence. 
This  case  is  actually  simpler  than  that  discussed  here  and 
is  analyzed  in  Appendix  F  of  Ref.  19.  The  second  case  is 


on*  pair  of  complex  conjugate  *ig*nvalu*».  This  cam*  is  of 
little  practical  interest  here  and  to  this  author's 
knowledge  has  never  been  analyzed  from  an  analytical  point 
of  view. 


g..£  Mail Mn MT  Flutter  Solution 


For  the  sake  of  simplicity.  w*  will  assume  that  the 
second-order  nonlinear  terms  are  equal  to  zero  (this 
corresponds  for  instance  to  a  “symmetric11  system,  i.  e.  .  a 
structure  that  is  invariant  to  changes  in  sign  of  the  state 
variables).  As  shown  in  Refs. 11  and  15.  the  case  of  nonzero 
second-order  nonlinear  terms  adds  complexity  to  the  analysis 
without  adding  any  new  feature  to  the  phenomenon.  An 
approximate  solution  of  Eq.  i  can  be  obtained  as  follows.  Let 


so  that 


with 


A  *  Ae  4  a*  \a 

\t  7  ±  I 


(10) 


(11) 


(12) 


and  assume  that  the  solution  for  x  is 
follows: 


where 
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T  * 


a  function  of  time  as 

(13) 

(14) 


This  means  that  time  dependence  appears  in  two  different 
forms:  the  first  on*  is  independent  of  fr,  while  the  second 
one.  the  so-called  slow-scal*.  becomes  slower  as  e  goes  to 
zero.  On  the  other  hand.  introducing  the  transformation 
given  by  Eq.  14.  the  dependence  of  x  upon  x  is  independent 
of  e.  whereas,  the  dependence  upon  "the  fast  scale.  t/€*. 
becomes  faster  and  faster  as  &  goes  to  zero.  Not*  that  t  and 
T  correspond  to  the  variables  t«  and  ta  in  Refs.  11.12.14 
and  13.  Because  we  are  limiting  our  analysis  to  0(fr*).  only 
these  two  time  scales  are  needed.  In  this  sense,  the  method 
used  her*  (multiple  time  scales)  is  different  from  the  two- 
time-scale  method  used,  for  instance,  by  Col*  and  Kevorkian 
(Ref.  13). 


Equation  13  implies  that 


<L  = 

dJc 


2 


+  * 


dr 


W*  will  also  assufflo  that  A  and  c 
of  A  so  that  ~  "IT 


(15) 

are  analytic  functions 


A(\)  s  A0  +  LX  XA  Aa  +  0(t4) 

whore  A0  ■  A<  >^)  and  Az  -9^3 ^.whereas 

)  ”  C*Ttf  r‘^  A  #  )  f  0(fi  ) 
Consider  a  solution  of  the  type 


(16) 

(17) 


X  s  £  X,  •+  fi3  Xa  +  O 


(IB) 


where  x^  are  functions  of  t  and  T. 

Combining  Eqs.  1  and  15-18*  yields 

l  [#*./»*  -  6.  !,] 

+  t*  [«*,/«♦  +  f»] 

+  o  c*r)  =  o 


where  all  the  terms  in  brackets  are  not  explicitly  functions 
of  £  (it  may  be  noted  that  this  is  the  mathematical  reason 
for  introducing  the  slow  time  scale  |)  and 

U  -  f  *Z  c»ffr( *')  K'f  *t-  }  (20) 

Separating  terms  of  the  same  order  of  magnitude  yields 


3*,/H  -  A0  x,  s  o 

3x4/at  -  Ao  *i  ^  -»  4  f* 


Disregarding  the  terms  in  the  solution  corresponding  to 
characteristic  roots  with  real  parts*  the  solution  of  Eq.  22 

4*  M 

x,  =  a(x)  u  c  +  <C(x)  a  t  (23) 

where  i<^  is  the  imaginary  eigenvalue  of  A.  and  u  is  the 
corresponding  eigenvector  (see  Eq.  B):  the  main  difference 


between  the  solution  of  tho  linear  problem/  Eq.  7*  and  the 
first  order  solution  of  the  nonlinear  problem*  Eq.  S3*  is 
that  c  is  a  constant  whereas*  in  general*  a  is  a  function  of 
T  (note  that  in  Eq.  SI  the  time  derivative  is  a  partial 
derivative*  whereas  in  the  linear  problem*  Eq.  3*  the  time 
derivative  is  an  ordinary  derivative). 


Substituting  Eq.  S3  into  Eq.  S3  and  noting  that 

f ,  =  c  M  ^  -fl'^+CCT) 

t  ft' 

(a  y  «,  T) 

-fr  +  CCT 


where  CCT  indicates  complex  conjugate  terms  and 

C-  £  c Vf  I  <«> 

Ur  j 

C*pjr  U*"*  UfU1Ur)  (36) 


one  obtains 


O  /o  L  1  1  J  ^  ^*4 i 

+  u  9o  /0  t  4  m  a  +  -f j  «V]  e 


•COjfc 


+  CCT 


The  solution  contains  "secular  terms"  (i.e. *  terms  of 
the  type  teu**  )  unless 

Ci )  i 

axt  w  9a  /9r  +  *4 ,,  u  a  4  ^  aV -  0  <28) 


_T 

where  -  is  the  left  eigenvector  of  A,  corresponding  to  the 
eigenvalue  i <  i.  e.  >  the  nontrivial  solution  of 

~T  fiw.  X-  A.  ]  *  0 

L  "  (29) 

If  v  is  normalized  so  that 


AT T  U  =.  { 

Eq.  28  may  be  written  as 

+  £  «■  +  %  s  0 

c  ' 


3 


where 


?--  -  X2  ?T  Ai  U 

-  -i  f3 


By  setting 


A  *  |*| 

and  separating  real  and  imaginary  parts  one  obtains 


(30) 

(31) 

(32) 

(33) 

(34) 


d\a\  /0x  +  |*|  +  hi3*  0 

v<f/0 t  +  jSj  +  yr  |«T  »  0 

Equation  39  may  be  easily  solved  by  setting  lal 
to  obtain 

0 


or 


.  -  2L(i  +  ) 

*>«• 

(where  k  depends  upon  the  initial  condition)  and  hence 

l/i 


A  |  ' 


-  |g«  /  y* 
UM  •2/s«-r 


(35) 

(36) 
-i/VT 

(37) 

(38) 

(39) 
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g 

t 


I 


Substituting  into  Eq.  36  and  integrating  yields 


i  —{( ft  *  ^  mv  dx  *  i. 

o 

=■  fai ?«./u-?*)x  + fa /?*)■(«  I*  I  *  <f. 


(40> 

Combining  Eqs.  10.  14.  18.  23.34.39  end  40  one  obtains 


X  |a|  «)*0(t,)4it 


with 


W  r  <J#  +  l  +  J  A  -  j 


(42) 


(vlt.r.  H= 

f.  '(/./?.  )-Ma  I  ■»  V. 


(43) 


whereas  iu)#  is  the  imaginary  eigenvalue  (and  y  the 
corresponding  eigenvector)  of  A0  (see  Eqs.  8  and  9)  and  lal 
is  given  by  Eq.  39. 


2J1  ttimtl.fl.flP  ftf.  Solution 


*(x 


The  function  lal  divided  by  |^/^|  *  i- »  (see  Eq.  39) 

I*  I  »  s^»v  (p*  /y«.) _ 

r.) 

(44) 


I  I 


•/l 


where 


I-  i ^ 


h  e 


T.  =  -  I  k  |  /«  /5ft 


(45) 


is  plotted  in  Figs.  1-4  for  all  of  the  possible  combinations 
of  the  signs  of  and  Only  real  values  of  lal  are 
plotted.  It  is  worth  noting  that  the  absolute  value  of  k 
influences  only  the  value  of  x,  which  in  turn  determines  the 
value  of  the  abscissa  that  corresponds  to  the  time  origin  t » 
0.  Note  that  (see  Eq.  32) 


f «  »  0  *’r  A.:11  (A>‘M 


(46) 
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since  the  linear  system  it  unstable  (i.e.  .  /gft< 0)  for  1 

(i.  o.  .  A  >  A,  ).  Consider.  separately  tha  cases  ^>0  and  ^<0. 

£±11— I— Ljk2fll 

Considtr  first  Fig.  1.  which  corrtsponds  to  >  0 
(stabilizing  nonlinear  terms)  and  A  >  X0  (i.e.  .  0# 

destabilizing  linear  terms).  In  this  case*  there  exist  two 
branches!  depending  upon  the  sign  of  k.  If  the  initial 
conditions  are  such  that  at  time  t  "0/  lal  >  'then  k 

<  0  and  the  amplitude  decreases  to  the  limit-cycle  value 

^  M  s  lalLC  s  (47) 

T  -*•* 

Viceversa  if  lal  >  I  /&*./  Vr  l'**  #  then  k  >  0  and  the  amplitude 
increases  to  the  same  limit-cycle  value.  In  other  words#  for 
yR>  O  and  X>Xa  (i.e.#  fcC  0)#  there  exists  a  stable  limit 
cycle;  note  that  the  solution  lal*0  is  unstable.  This 
behavior  could  have  been  predicted  from  Eq.  31  by  noting 
that  for  large  values  of  lal  the  stabilizing  nonlinear  terms 
are  dominant.  whereas  for  small  values  of  lal  the 
destabilizing  linear  terms  are  dominant. 


Next  consider  Fig.  2#  which  corresponds  to  yR>  0 
(stabilizing  nonlinear  terms)  and  A  <  A,  (i.e.#  J&R>  0# 
stabilizing  linear  terms).  In  this  case#  real  values  for  lal 
are  obtained  only  if  k  <  O  and  T .  The  solution  tends  to 

l*ro:  #,  i  i  « 

JLr*.  |  «.  |  *  0 

t-*—  (48) 
i.e.#  the  solution  lal  -  0  is  stable.  This  behavior  is  to  be 
expected#  since  both  linear  and  nonlinear  terms  in  Eq.  31 
are  stabilizing. 


Summarizing,  if  7R  >  0  (stabilizing  nonlinear  terms), 
there  exists  a  stable  limit  cycle  for  A  >A©  (when  the  linear 
analysis  predicts  instability)#  whereas  for  X<  A0  the 
solution  lal  *  0  is  unconditionally  stable. 


C«»f  S  L&SQ2 

Next  consider  Fig.  3#  which  corresponds  to  ?R<  O 
(destabilizing  nonlinear  terms)  and  X  >  X0  (i.e.#  O. 


destabilising  linear  terms).  In  this  case*  lal  is  real  only 
for  k  <  0  and  T  <  Te  .  The  solution  tends  to  infinity  as 
-t  tends  to  r  : 


i.  e.  •  the 
behavior  is 


I  A  I 

T  ^  TV 

solution  is 
predictable 


I  a.  I  *  oe 

1  (49) 

is  unconditionally  unstable.  This 
ble  since  both  linear  and  nonlinear 


terms  in  Eq.  31  are  destabilizing. 


Last  consider  Fig.  4>  which 
(destabilizing  nonlinear  terms)  and 
stabilizing  linear  terms).  Note  that 
k  %  0,  respectively.  In  both  cases 


corresponds  to  fa  <  0 
A  <  A,  (i.  e.  ,  pn>  O. 
two  branches  exist!  for 


I  al  •  =  Ipt/frl 


However<  if 
at  time  X 
infinity: 


x 

the  initial 
-  0  (i.e. 


conditions  are 
k<0)i  then 


such  that  lal  >  la 
the  solution  grows 


T-*T, 

On  the  other  hand»  if 
lal  <  laluc  at  time 
goes  to  zero: 


I  A  I  =  0  '  ^  la  I «  ^  |  ^  I LC 

1  (92) 
The  preceeding  behavior  may  be  character! zed  by  saying  that 
there  exists  an  unstable  limit  cycle.  The  solution  lal  *  0 
is  stable*  whereas  the  solution  lal  ■  lalLC  is  unstable.  In 
other  wordsi  the  system  is  conditionally  stable  (i.e.!  is 
stable  only  if  the  initial  conditions  are  sufficiently 
small!  whereas  it  is  unstable  if  the  initial  conditions  are 
sufficiently  large).  This  behavior  could  have  been  predicted 
from  Eq.  31  by  noting  that  the  stabilizing  linear  terms 


^  /*/, 

(91) 

initial  conditions  are  such  that 
0  (i.e.!  k>0)!  then  the  solution 

0 


(92) 

characterized  by  saying  that 
cycle.  The  solution  lal  *  0 
lal  ■  lal..  is  unstable.  In 


dominate  for  small  values  of  lal,  wherea 
nonlinear  terms  dominate  for  large  values  of 
be  emphasized  that  this  behavior  is  very 
practical  applications!  since!  if  the  initial 
sufficiently  large!  instability  occurs  even 
analysis  predicts  stability.  The  primary 
flutter  taming  is  to  prevent  this  phenomenon. 


whereas  destabilizing 
ues  of  lal.  It  should 
very  significant  in 
nitial  conditions  are 
even  if  the  linear 
rimary  objective  of 


agaaflaE 


Summarizing*  if  TfK<  0  (destabilizing  nonlinear  terms)* 
there  exists  an  unstable  limit  cycle  for  A  <  A* (when  the 
linear  analysis  predicts  stability)*  whereas  for  A  >A0  the 
solution  is  unconditionally  unstable. 


In  conclusion*  the  results  of  this 
restated  as  follows:  there  always  exists  a 
by  (see  Eq.  41) 


section  may  be 
limit  cycle  given 


j  =  2  |U-x.)  eK  /y„  | (a 

1  r  ...  (»3) 


where  «/>*  -  (Yt  *•  •  constant.  If  **>  0 
(stabilizing  nonlinear  terms)*  a  stable  1 imit  cycle  exists 
for  A  >A*  (with  amplitude  growing  like  > »  whereas  the 
system  is  stable  for  X<  Ao  •  This  behavior  is  indicated  in 
Fig.  3*  where  the  arrows  indicate  the  variation  with  time. 
Viceversa  if  #*<  0  (destabilizing  nonlinear  term)*  an 
unstable  limit  cycle  exists  for  A  <  A*  (with  amplitude 
growing  like  /vX)  whereas  the  system  is  unstable  for  A>A.. 
This  behavior  is  indicated  in  Fig.  6. 
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SECTION  III 


THEORY  OF  FLUTTER  TAMINO 


UL  In&mAiuLtlan 


The  results  of  the  analysis  of  nonlinear  flutter 
presented  in  Section  II  indicate  -that  the  behavior  of  the 
aeroelastic  system  beyond  the  flutter  boundary  (e. g. *  benign 
versus  destructive  flutter)  is  determined  by  the  sign  of  the 
real  part  of  a  complex  number  ?  defined  by  Eq.  33.  This 
suggests  a  simple  procedure  for  flutter  taming:  modify  the 
system  so  that  the  additional  contributions  produce  the 
desired  positive  value  for  fa  .  It  is  convenient  to  add  the 
requirement  that  the  linear  characteristics  of  the  modified 
system  be  equal  to  those  of  the  original  system.  This  is 
important  in  order  to  avoid*  for  instance*  that  the 
introduction  of  flutter  taming  decrease  the  flutter  speed. 


It  may  be  worth  recalling  that  the  objective  of  flutter 
suppression  by  active  control  is  to  modify  the  linear  system 
to  achieve  an  increase  in  the  flutter  speed  (Ref.  3)  the 
above  requirement  that  the  linear  characteristics  before  and 
after  the  flutter-taming  modifications  ensures  that  flutter 
taming  may  be  made  on  a  system  already  modified  for  flutter 
suppression  without  affecting  the  flutter-suppression 
modifications. 


Two  types  of  modifications  are  considered  here.  In  the 
first  one  (type  A).  the  nonlinear  terms  of  the  system  are 
modified  to  ensure  the  benign— flutter  condition.  These 
modifications  could  be  a  modification  of  the  nonlinear 
elastic  characteristics  or  the  addition  of  a  nonlinear 
damper.  In  the  second  one  (type  B>«  new  state  variables 
(typically  active  control)  with  a  nonlinear  feedback  are 
added  to  the  system  (augmented  system).  This  can  be 
accomplished*  for  instance*  by  having  a  sensor  (e. g. *  an 
accelerator)  and  an  actuator  which  controls  a  control 
surface*  with  a  driving  force  proportional  to  the  cube  of  a 
linear  combination  of  the  original  state  variables. 


Whereas  many  different  modifications  may  be  included* 
only  the  above  two  types  are  considered  here.  The  modified 
system  is  then  described  by 


04) 


dx/dt  =  AM  x  +  *)+*'/&,*) +  BU)y 

dy/d*  D  y  +  t( cTx)J 

The  term  *'f'(x»A)  (where  is  a  constant  and  f'  it  a 

nonlinaar  vactor)  corresponds  to  the  type-A  modification.  On 
tha  other  hand*  the  term  Bu  (with  governed  by  Eq.  55) 

corresponds  to  type— B  modification.  The  nonlinear  feedback 
is  provided  by  the  term  kmz(ctx)3  (where  z  and  c  are 

arbitrary  vectors  and  is  an  arbitrary  constant).  Note 

that  (as  shown  in  details  later)*  the  lack  of  linear 
feedback  in  Eq.  55  guarantees  that  the  linear 

characteristics  of  the  modified  system  are  equal  to  those  of 
the  original  system.  Note  also  that  by  setting  *'  ■  O  one 
obtains  the  modification  of  type  B *  whereas  setting  B  ■  Q 
one  recovers  the  modification  of  type  A.  For  the  sake  of 
convenience*  the  analysis  is  carried  out  under  the 
assumption  that  both  modifications  are  present  at  the  same 
time. 


SLR  ElilUlX.  lamina  Solution 

Equations  54  and  55  may  be  recast  as 

dx+/cU  *  A*(M  **  +  A) 

where 


(56) 


(57) 


Equation  96  is  formally  equal  to  Eq.  1  and  therefore 
the  results  of  Section  ZX  are  applicable  to  Eq.  96  as  well. 
Hence*  the  solution  of  the  linear  system  is  obtained  by 
setting 

r-  **]  * 

*  *u[‘!-V]°o  <60) 

which  shows  that  the  eigenvalues  of  the  modified  system  are 
equal  to  those  of  the  original  system  plus  those  of  the 
actuator  system. 


si- A  -B 

9  s  X-p 


In  particular*  for  an  eigenvalue*  s  -  icJ  »  is  on 
the  Imaginary  axis  and  (assuming  that  the  actuator  system 
does  not  have  an  eigenvalue  equal  to  i ui#>  the  corresponding 
eigenvector  m+ *  i. e. •  a  nontrivial  solution  of 


i -a;]  «+ 

-B, 

O  tloJ.I-D  J 

(where  A#  ■  A(  *,)  and  B#  »  B(  A#>  >  is  given  by 


uf  =  O 


u 


+ 


(61) 


(62) 


(where  u  is  the  eigenvector  of  A. 
addition*  the  left  eigenvector*  aT*» 


*  defined  by  Eq.  9).  In 
i.  e.  *  the  solution  of 


A0  -  B# 

I  ~  P 


is  given  by 


where  v  is  the  left  eigenvector  of  A# 


and 


(63) 


(64) 


aTt  s  yT  B  I  -  .J)  J 

The  condition  for  benign  flutter  is  given  by 
0*  *  >  0 


(65) 


(66) 


with 


r*rT  n 


<o 


(67) 


where 


with  f 


o) 


,(•>  +i  +  *  ‘a 

fa  s  .  , 

,  K"  *  3 J4 yU# 


(68) 


,  given  by  Eq.  26  end  f'j *  given  by  a  similar 

T 


expression*  whereas 

/•*  5'  !f 

Using  Equations  64  and  6S>  Eq.  67  may  be  written  as 

t'  -  y  +  +  «•  y" 

where 


(69) 


(70) 


y  *  *T  is 


<•) 


(71) 


is  due  to  the  nonlinearities  of  the  original  system 


l  t  f 1  ^ 

y  *  ?  f» 


(72) 


is  due  to  the  additional  nonlinearities  of  the  control 
surface  aerodynamics  and 


r«VXiT* 

The  condition  for  benign  flutter  is 

u  *  ar*  +  K,3ri  + 


(73) 


(74) 


which  indicates  that  if  jr*  1*  0  (or  ^  0).  it  is  always 
possible  to  choose  %•  (or  *”)  to  satisfy  the  above 
condition.  Note  that  it  is  easy  to  ensure  that  f  O.  For* 
it  is  always  possible  to  choose  c  such  that 


|>^|«  |&T  *|  $  0 


(75) 


and  to  choose  z  such  that 


<  -  \r\‘  (r$T*)  +  °  (74) 

As  a  final  remark!  not  that  as  t«'sgn(  )  goas  to 

infinitg<  goas  to  infinity  and  hanca  tha  limit-cycle 

amplituda  goas  to  zero  (similarly  for  *">.  In  othar  words! 
tha  amplituda  of  tha  limit-cycla  can  ba  mada  as  small  as 
dasirad  by  choosing  appropriata  valuas  for  tha  constants 
and 


SECTION  IV 


AIRFOIL  IN  SUPERSONIC  FLOW 


4.  1  Qovrnlno  Equation* 


Consider  an  airfoil  in  a  supersonic  stream  (see  Fig. 
B.  1).  Let  the  airfoil  have  only  two  degrees  of  freedom: the 
vertical  displacent* h  ■  b£ (where  b  is  the  semichord >* of  the 
center  of  mass  and  the  rotation*  oc  ,  around  the  center  of 
mass.  Let  the  airfoil  be  restrained  by  two  springs  at  a 
point  E  (representing  the  elastic  axis  of  a  wing).  The  first 
spring  produces  a  vertical  force,  function  of  the  vertical 
displacement  at  E.  The  second  produces  a  moment  function  of 
U  .  Let  the  two  springs  be  of  symmetric  nature  (e. g. *  the 
same  force  is  obtained  in  tension  and  compression)  and 


defined  by  the  constants  KB  « K*  » K,'  *  and  K*  *  see  Eg.  B.  41. 


Assume  that  the  piston  theory  gives  a  sufficiently  accurate 
description  for  the  pressure  distribution.  Then* as  shown  in 
Appendix  B*  the  governing  equation  are  given  by  Eq.  B.  91 


(77) 


A  = 


1  be 

recast 

as  Eq. 

1  with 

XT 
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0 

0 

i 

0 

0 

0 

0 

0 

0 

0 

0 

i 

0 

0 

0 

0 

*  fc’ll 

-drS« 

9S, 

0 

5. 

0 

0 

~ku 

-  ICn 

«s, 

0 

-•Si_ 

0 

0 

V 

V 

V* 

V 

and 


(78) 


4  ; 


0 

0 


(A) 


Cs) 

J*-L 


,(a) 

* 


(79) 


In  Eq.  78  v  is  the  nondimensional  moment  of  inertia  (see 
Eq.  B.  19),  k-  are  constants  (given'  by  Eqs.  B.  48  and  B.  44  in 
terms  of  the  spring  constants*  the  speed  of  sound  aM  and  the 
mass  m)*  8^  are  constants  given  by  the  Eq.  B.  33  in  the  terms 
of  (nondimensional  location  of  midchord  with  respect  to 


;  j-.  -r. 


the  cantar  of  mac*).  In  addition  C  it  th*  matt  ratio  (givan 
bg  Eq.  B.  27)  and 

A  s  Md  a  ii»  _^*P««4  (BO) 

(whara  M»U,/o.,,  is  tha  Mach  numb ar)  is  proportional  to  tha 
spaad  of  flight. 

In  addition  tha  structural  nonlinaar  tarms  ara  givan  by 
(saa  Eqs.  B.  44.  B.  49  and  B.  30) 

(*)  =  j_  k.S«s-  k'  +  f*)s 

3$  6  *  5 

j*’  (s )  =  -  H  *  $ «' V§ k^-  O’1*-  **)’ 

(81) 

Finally  tha  aarodynamic  nonlinaar  tarms  ara  givan  by 
(saa  Eqs.  B.  33»  B.  36«  B.  37»  B.  39  and  B.  40 

<j%)  *  «  (c„  s,  +  ca  S  «*4  4  =„  S, 

+«x(s.v,J  +  3S'v'li  +  3!vJ*a  + 5siJ) 

{?(*-)  * «  (**  s,  n  «» 4  q  -*  i  +  c„  4 
+6£  (q  +  34  ■'•J1*  +  34*0**’-* 4**^  <ez> 

with 

-UTs  M  0/  -  $  (83) 

wharaas  (saa  Eqs.  B.  33  and  B.  36) 

r.*/3  1  -'/*■' 

-,/z  °J 

Thasa  valuas  ara  obtainad  by  assuming  that  tha  downwash  usad 
in  tha  piston  thaory  is  in  tha  diraction  of  tha  rj-axis. 
Similar  txprassions ara  obtainad  by  Smith  and  Morino  (Rtf. 
IS)  by  assuming  tha  downwash  to  ba  in  tha  diraction  of  tha 
y-axis:  in  this  cast  tha  coafficiants  ara  givan  by  ittt  Eqs. 
B. 39  and  B.  40) 
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(89) 


,  ■ 


i/tf  -i/a 


;] 


The  solution  is  of  tho  typo  given  by  Eq.  9>  where  s{  and 
U;  aro  oigonvaluos  and  tho  corresponding  eigenvectors  of  A. 
In  tho  caso  under  consideration  here. the  oigonvaluos  aro  tho 
solution  of 


AS.  _*S, 


Ki«.  ktt+\&± 

IT  v 


s  tA 


or. adding  third  and  fourth  columns  (multiplied  by  4  )  to 
first  and  second  columns  rospsetivoly  and  multipling  tho 
fourth  row  by  v» 


Sitt 


s*  +  *  S0  a  +  k„ 
.si  s+k. 


-  o  S,  a  +  k-,x  -  X  S0 

S1-  0^54-  4-  A 


which  yiolds 


4^  5*  4-  A4  S J  S  *  +  4,  S  f  A.  *  0 


with 


a4*  V 


0-(V  S.  4  S4) 

«te  V»  K||  +  ku  4  ) 

+  Sj.  If,,  .2  5,  fen)  (89) 

Ifii  few  —  kr*  4  A(s,  K%  i  ■+•  S„  fci*) 

At  tho  stability  boundary.  X  •  A,,  a  pair  of  complex 
conjugate  roots  aro  puroly  Imaginary  (s  ■  t  i«*»# .  see  Eq. 
4). and.  separating  real  part  and  imaginary  part. one  obtains 

a.f  UJ*  -  4^  oJ*  +  4o=  O 


^4  <0#  -  s  ® 


or.  from  Eq.  91. 


(92) 


,\  .VkV  V^  .MV-V.n  *-  •-  -  * 


<0/  -  (5ft  trlfc  +  4  k,t  -  2St  k,*) /(»$*+  **) 

and  1  from  Eq.  90# 

A  VO).4-  (W  +  fr,  Kx  -  ♦ «' VcwgJ 


-u^4  *S(  +  S(  fc>i  •*  S„  #f, 


a 


C93) 


It  nay  ba  worth  noting  that# by  using  tho  expression  for  k^r 
and  (Eqs.  B.  48  and  B.  33)#  Eqs.  92  and  93  may  bt  rewritten 

a« 


<94  > 


and 


y  .  .  CO/-  Kt)«  -  >M*)  ♦  o/  >  cV/3 

“>.*  €.  ♦  (S-i.)Ks 


(95) 


in  agreement  with  Ref.  18  and  Ref.  15  (where  kt  mfl$  <  k«* 

£,«  «7t  and  <T-£«»  ■  *?■  mote  a  misprint  in  Eq.C16  of  Ref.  15 
where  V  should  be  <o/). 

The  eigenvector  u  (the  nontrivial  solution  of  Eq.  8)# 
normalized  by  setting  ^  ■  1# is  given  by 

i  % 


U  - 


u. 


(96) 


where 


M.  =  (r  4  1  tJ#  ®  S0  4-  ku  «  £,  +  k,a  -  A,  &,  ) 

e  &  £t  +  k,x  )/ (-  V  +  “J0  6  ^  ) 


(97) 


Similarly  the  left  eigenvector  #\f  is  given  by 

aTTs  J-  L*0.  ^ 

*0 

where 


(98) 
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s-s* 


K  *  «d 


(99) 


and 


aTs  (.  <j/  4i^cS(  +  k»i)/ ( + 

*  (-  lu>»  °^i  +  “  \^o  )/(-  *  “*4  c  +  ) 


/at,  J  Uh,  +«4  -»*,  *1  r  I  ) 

hM  [  ,u0vf«^J  {^) 


(100) 


whereas  o(,  it  a  normal isat ion  constant  such  that  <irTU  >  1  (saa 
Eq.  30).  1.  a.  . 


*#z  2!oJt  (I+U'V)  +  «•  [4  +fr+'')3  *  uv-44]  (101) 


4.3  Nonlinaar  Flutter 


Using  Eqs.  20  and  24  ona  obtains* for  tha  spacific  casa 
undar  consideration.  , 


with 
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?53a 

-  V6  k.  <T(*,u.  etM*fc 
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-  ^  pL(,|  +  $M.Je'uifc  +  cctJ3 
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+  Cft  (4  M  <  ’^+  CCT )  (ftii*;,  «  *4-  cct) 

+  s'  c,s  ccT)a(au»ij# c1*0***  cc r) 

4 S„  e'^r  ccr)3 

4<r^3S,  (  *^fc  +  cct)2(’a«  i’«,e'Ht4  ccr) 

3Si  (***  C.  CCT)(*U  .’cj,  e<U**+  cc  t)2 
je-us.  LUJ  < M0#*+CCr)1 


(103) 
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(*.  ccr  )(auc,u,\  cot)* 
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+  *  ^  £.  ^^e^^ccT)1 

4  ^  s,  (a  U.  C,UJ#t+  CCT)1  (a.  »'ccJ„  c“^fc+  ccr) 

+  r  CM  5X  (ftwe'^^ccT^AKt'^^Vccr) 

**XS/  e'U*6  +  CCT-)J 

4*X  3Sa  C**0  *  "d‘t+  CCT)a(au  iw,  €  ;<**+CC r) 

4  ^3Sj  (o.^  «. ,ui,t  +  ccr)  (auw0  €  ,l**>t  +  ccr)z 


where* CCT  indicate*  complex  conjugate  term*  and 

'Cr  s  Mat  i(0( 

Equation*  103  and  104  may  be  written  as 
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i  # 

r») 

*  *  e 
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where 

?Ci  *  ^  *  3  u*“*  3  (U  S*)1  (!+&«*) 

+  «■  [c„  4  n*'*l«*+  c(X  s,  (3uH^-Ma;tio.+clS  s,  «V 

T®x  [S.-3^**3*  +35,  (Swv*  u.  J  tij0 

+  3&a  (to  2mm*  *o  *  w4)  u^2,  +<4^  3w.aw*iV^J 
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f'l  =  «/*  ktifruu** 

-  ^  3(t*Ju)*(l  +  Suy 

4  ^  [cai  &,  H  3  u*u*+ cxt  St  (a  u  u*~  2)  i<ot + cAi  $A  u*«*  iwc 

+*%  [$,  3*0*3*  +  3^(i?w*51ru  -  +jzu*)  iK)0 

+  3  (vJ  &  *  w*  -  3  M1)  Mi1,  +  So  3  u*m*  /u>,3J  ( l08) 

Similar  expression  for  g  '  and  g^  art  not  explicitly  given 
here  because  they  are  nos  needed  in  the  current  analysis. 

The  condition  to  avoid  secular  terms*  Eq.  29.  yields  Eq.  31 
where  now  (Eqs.  78* 96  and  98) 

/  f  a  r\  a  ^  T  /  J  "s 


P  s  -  \  a  s 

v  a/ 
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xtO  U 


-.it  I  I  vvl  r  Lai(5.-^S)h/< 

5;  o  -s,/v  l*-j  1  1  '  “<>’> 

and 

{,»  ( “W  +v^iA* 

•  /  *  )  [  #5»  \ 

^  jJ*/V  '  '  <“<>> 

The  limit  cycle  amplitude  is  then  given  by  (see  Eqs.  47  and 
90) 

I  e«4f  $0U-  AT  S.  u  )/*.!  I 


^  [(^i  f  v  )/*•] 


The  numerical  result*  ere  discussed  in  Section  V. 


..Clttttet-JffW  till 


In  order  to  implement  the  f lutter-taming  formulation* the 
system  considered  above  must  be  modified.  Three 
modifications  are  discussed  here  (the  details  are  given  in 
Appendix  B)  The  simplest  one  consist  in  changing  the 
nonlinear-stiffness  constant  k*  (this  case  “simulates"  a 
modification  of  the  structural  design  of  the  aircraft  to 
insure  benign  flutter).  Another  simple  approach  consists  In 
adding  a  nonlinear  damper  (this  case  "simulates"  a 
modification  of  existing  aircraft).  These  two  approaches  may 
be  studied  using  the  same  equations  used  above  where  now  g*. 
is  given  by  Eq.  61«or 

h  *'  **  <>»=> 

This  implies  that  V*is  given  by  Eq.  70  where  O  and 
(see  Eq.  72)  0 

s  -  IcJ*  3  UXU*  (113) 


In  a  more  sophisticated  approach  one  can  make  use  of  a 
control  surface  and  nonlinear  active  control.  The  presence 
of  the  control  surface  Introduces  forces  and  moments  (where 
0  is  the  deflection  of  the  control  surface).  The  equation  of 
the  actuator  dynamic  is  assumed  to  have  an  input 
proportional  for  Instance  to  £J  (e.  g.  #  the  doubly-integrated 
signal  from  an  accelerometer  located  in  C). In  this  case  the 
governing  equations  are  given  by  Eq.  B.  69  which  may  be 
recast  in  the  form  of  Eq.  94  with  4c’*  0  (see  however  Eq.  112) 
and 
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This  yitldt  (see  Eq.  65) 
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.  f-  S*  +  x  (-$<**!)]/(&  + 


120) 


V 

and  a  similar  expression  for  v4  which  is  inessential  here. 

The  expression  for  $+  is  given  by  Eq.  70  where  ?'  is  given 
by  Eq.  113  (with  x’  if  the  damper  is  not  added)  and  j^given  by 
Eq.  73«  or 


3*"  '  3/*V*  5T  *  s 

v 

with  *T  given  by  Eq.  120  and  given 

/•*  ST*  -  A 


v>,  * 

by  Eq.  69  or 


(121) 

(122) 


combining  Eqs.  120.121  and  122  one  finally  obtains 

a "  *  3  [^.r(-^‘"+  VS")+  U23) 

In  the  design  of  the  actuator. the  parameters  J.r  end  k  may 
be  used  to  insure  that  0  (or  even  2^x.s  optimal) 

Numerical  results  are  considered  in  Section  V. 
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SECTION  V 


NUMERICAL.  RESULTS 


9. 1  Introduction 

Jn  order  to  assess  the  feasibility  of  flutter  taming 
(e. g.  ,  whether  the  modifications  may  be  implemented  with 
physically  acceptable  values  for  the  constants  -k'  and  ■k*  )« 
specific  numerical  results  are  presented  here.  All  the 
results  are  obtained  using  a  configuration  similar  to  that 
of  Ref.  19  and  18.  In  particular  the  following  values  are 
used  for  the  constants  ,  S  ,  v  ,  fc5  and  K's 

€o-“-2 

-  ££- .2 
v  -  .  21 
ks-  0. 

kE«  0.  .  (124) 

9.2  Flutter  Analusls  of  Original  Sustem 

Consider  first  the  linear  system  according  to  the 
analysis  of  Section  IV  the  •  flutter  frequency  is  given  by 
Eq. 94  or»  using  the  above  values  for  K and  §  tone  obtains 

5  * 


l 

.  21  4  1/3  4.  04 


« 


U 

7 


.  36  Jl* 


(129) 


where  (as  in  Ref.  19) 

a.  *  a*/* 


(126) 


The  flutter  boundary  is  given  by  Eq.  99  and  is  plotted 
as  a  function  of  Vft£ln  Figure  7  (for  1/e  ■  4#  St  12  and  16). 
The  corresponding  Mach  numbert  M0  ■  A*/*  (see  Eq.  80)t  is 
plotted  in  Figure  8  as  a  function  of  *££(for  the  same  values 
of  9  ).  The  above  results  coincide,  as  expected,  with  those 
of  Ref.  19. 

In  the  nonlinear  analysis  two  formulations  have  been 
presented.  The  two  differ  only  for  the  direction  of  the 
downwash  used  in  the  piston  theory.  In  the  first  onct  the 


downwash  is  in  the  vertical  direction  (akin  to  the 
formulation  of  Ref.  13):  1  will  refer  to  this  formulation  as 
type  V  (for  vertical).  Xn  the  second  one. the  downwash  is  in 
the  local  direction  of  the  normal  to  the  airfoil  surface:  X 
will  refer  to  this  formulation  as  type  N  (for  normal). 
Equation  83  is  used  for  type  V.  whereas  Eq.  84  is  used  for 
type  N.  The  value  of  the  limit-cycle  amplitude  (as  given  by 
Eq.  Ill),  are  plotted,  in  Figures  9-12  (for  1/6  ■  4.  8.  12 
and  16  respectively),  as  a  function  of  for  several  values 
of  ^(values  of  in  parenthesis  indicate  an  unstable  limit 
cycle).  The  results  obtained  using  the  first  formulation 
(type  V).  shown  on  the  upper  side  of  the  figures  coincide, 
as  expected,  with  those  of  Ref.  13.  Those  obtained  using  the 
second  formulation  (type  N)  are  shown  on  the  bottom  of  the 
figures  and  indicate  that  the  difference  between  the  two 
approaches,  while  not  dramatic  is  still  significant.  From  an 
physically  intuitive  point  of  view  the  second  approach  is 
obviously  more  satisfactory  than  the  first  one.  However  X  am 
not  aware  of  any  rigorous  proof.  that  indeed  the  second 
approach  is  the  correct  one.  Therefore  in  the  following  all 
the  figures  include  both,  type  V  (top)  and  type  N  (bottom) 
results. 


Figures  9-12  indicate  that  as  k^  increases,  the  limit 
cycles  tend  to  become  stable.  This  suggested  the  first, 
simple  approach  to  the  flutter  taming:  modify  the  value  of 
the  nonlinear  stiffness  constant  kit  in  order  to  achieve 
flutter  taming.  The  results  are  ■  shown  in  Figures  13-16  (for 
1/tr  m  4,  g,  12  and  16  respectively)  where  the  limit  cycle 
amplitude  is  plotted  versus  ki/J^f or  several  values  of/icjvalues 
of  in  parenthesis  indicate  an  unstable  cycle):  t  as 

mentioned  above  the  downwash  expression  of  Ref.  13  (type  V. 
Eq.  83)  is  used  for  the  top  graph,  whereas  the  formulation  of 
type  N  (Eq.  84)  is  used  for  the  bottom  graph.  These  figures 
are  a  different  way  of  presenting  the  same  data  shown  in 
Figures  9-12.  Xt  is  apparent  that  positive  values  of 
tend  to  stabilize  the  limit  cycle)  for  instance,  on  Figure 
13  (6m  1/4).  for  values  of  'Q/Kl  =  3  and  for  all  the  curves 
except  for  1.3.  the  ‘  limit  cycle  amplitudes  are 

relatively  small.  Note  that  if  lalUfca  i  .then  the  value  of 
the  solution  will  be  of  the  order  This  means,  for 

instance,  that  if  A  exceeds  A0  by  1%  (i.e. .  if  the  flight 
speed  exceeds  the  flutter  speed  by  1 X).  then  the  vertical 
displacent  will  be  equal  to  .lb  (3%  of  the  chord).  While 
this  value  is  relatively  small,  it  should  be  emphasized  that 
further  reduction  in  the  limit  cycle  amplitude  require  very 
high  values  of  k*  .hard  to  obtain  in  an  actual  aircraft 


dj*ign.  In  addition  f or  €T-  l/B  (Figure  14)  tho  curves  f or 
»  1.2S»  and  1.5  ara  not  stabilized  avan  for  3. 

For  ffm  1/12  and  1/16  (Figures  19  and  16)  only  the  curve  for 
JtT  .  5  is  stabilized. 

5.4  Taming  bu  Domnina 

Equally  simple  (at  least  from  a  mathematical  modeling 
viewpoint)  is  the  taming  by  addition  of  a  damper.  The 
results  (see  Eq.  70  with  ■k"-0  and  y'  given  by  Eq.  113)  are 
shown  in  Figures  17-20  (for  1/tf  *4*  8<  12  and  16 

respectively)  where  the  limit  cycle  amplitude  versus  k'  Is 
plotted  for  several  values  of  :  values  of  JiZ,  in 

parenthesis  Indicate  an  unstable  limit  cycle  (type  V 

formulation*  Eq.  83,  is  used  for  the  top  of  the  figure*  type 
N  formulation  Eq.  84*  is  used  for  the  bottom).  It  may  be 
noted  that  positive  values  of  the  constant  tc'  (i.e. , 
positive  damping)  have  stabilizing  effects  (whereas  this 
might  not  come  as  a  surprise*  this  point  deserves  further 
attention  and  points  to  the  need  for  the  development  of  a 
mathematical  formulation  which  would  identify  under  what 
conditions  nonlinear  damping  terms  are  stabilizing).  The 
same  comments  made  for  the  taming  by  stiffening:  for  1/4 

(Figure  17),  small  values  of  the  constant  k'  yield  a 
stabilization  and  than  a  dramatic  reduction  on  the  limit 

cycle  amplitude:  however  further  reductions  require  large 

values  for  the  constant.  For  8  -  1/8*  1/12,  1/16  (Figures 
18-20)  similar  trends  ara  observed*  but  stabilization  is  not 
achieved  far  the  higher  values  of  k  . 


The  expression  for  the  limit  cycle  amplitude  for  a 
system  modified  with  the  addition  of  active  control  is 
derived  in  Sect.  IV  (see  Eq..  70  with  fl'-0  and  «"  given  by 
Eq.  73).  The  constants  J,  r,  and  k*  determine  the  dynamics  of 
the  actuator*  Eq.  B.  67,  and* «  mentioned  at  the  end  of 
Sect.  IV*  they  may  be  used  to  optimize  the  limit  cycle  value. 
This  idea  of  "optimal  design"  -was  brought  about  by  the  first 
results  obtained  by  choosing  certain  values  for  J,  r*  and  k. 
The  diagrams  Indicated  very  little  influence  of  the  feedback 
even  for  K*m  20.  It  turned  out  that  those  values  yield  a 
value  of  y"  very  close  to  the  imaginary  axis  so  that 
(if  tfi-o^the  constant  &  has  no  effect  on  ftt»Eq.  70).  By 
changing  the  value  of  the  phase  angle  of  -  co*-j  +  iwr  +  k*  we 
found  that  a  value  for  close  to  tRe  maximum  is  obtained 

for  T/>i  -  .33  and  k/J  -  .  3.  In  view  of  the  fact  that  the 
nondimensional  moment  of  the  Inertia  for  the  airfoil  is 
v  — .  21,  it  seems  reasonable  to  chose  J  -  .1  and  hence  r  - 
.033  and  k  -  .03.  The  results  are  plotted  in  Figures  21- 


24  (for  1  8.  12*  16  respectively )*  which  show  .the  .limit 

cycle  amplitude  versus  tc"  for  several  values  of  ^(values  of  ^ 
in  parenthesis  indicate  unstable  limit  cycle).' Note"  that 
this  time  only  the  curve  1.5  (for  e«  1/8*  1/12.  and 

1/16)  is  not  stabilized  at  4c"  “-1.  Otherwise*  as  in  the  case 
of  the  taming  by  stiffening  and  by  damping*  small  values  oj'-e* 
yield  a  stabilization  of  the  unstable  limit  cycle  and'a 
dramatic  reduction  of  the  limit  cycle  amplitude.  Slower 
reductions  in  the  limit  cycle  amplitude  are  obtained  with 
further  increases  of  the  gain  -JCH  :  however*  high  values  for 
this  constant  may  be  more  easily  obtainable  in  this  case 
than  the  case  of  taming  by  stiffening  or  by  damping.  Further 
research  is  recommended  on  this  issue. 
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SECTION  VI 


CONCLUDXNQ  REMARKS 
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A  new  tool  for  the  aeroelastic ion*  flutter  taming*  has 
been  introduced.  The  objective  of  flutter  taming  is  to 
insure  that  flutter  is  benign  rather  than  destructive.  This 
objective  is  accomplished  by  ‘modifying  or  adding  nonlinear 
terms  to  the  governing  'equations*  or  by  augmenting  the 
system  by  adding  state  variables  with  a  nonlinear  feedback. 
The  general  theory  (based  on  an  asymptotic  solution  obtained 
by  using  the  multiple-time-scaling  method)  is  applied  to  the 
case  of  a  supersonic  airfoil.  Taming  by  changing  the 
nonlinear  elastic  terms  and  nonlinear  damping  terms  as  well 
as  by  adding  a  simple  feedback  system  (a  control  surface 
with  input  proportional  to  the  cube  of  the  vertical 
displacement  of  the  center  of  mass)  have  been  analyzed. 
Analytical  expressions  for  the  solution  have  been  derived 
and  used  to  obtain  numerical  results.  In  many  of  the 
specific  examples  considered  the  results  Indicate  that  it  is 
possible  to  achieve  stabilization  of  the  unstable  limit 
cycle  (which  means  that  destructive  flutter  is  replaced  with 
benign  flutter).  Also*  a  dramatic  reduction  in  the  limit 
cycle  amplitude  is  obtained  with  relatively  small  values  for 
the  multiplicative  constant  available  in  each  case.  Further 
Increases  of  these  constant*  however*  yield  only  marginal 
additional  gains.  As  decreases*  however*  stabilization  is 
not  achieved  (within  reasonable  values  of  the  multiple 
constants)  for  the  higher  values  of  .  The  most  promising 
results  are  obtained  using  active  control. 

fr.  8.  Ftcppmndgtlpnt. 

In  view  of  the  positive  results  obtained  so  far  and  in 
view  of  the  great  potential  for  practical  applications*  it 
is  recommended  that  this  work  be  continued  by  including  the 
influence  of  the  other  parameters  (such  as  £$*£  and  *  see 
Eq.  124).  It  is  also  recommended  that  the  work  be  continued 
with  a  wind  tunnel  simulation  to  confirm  experimentally  the 
results  obtained  here  (in  particular  on  the  feasibility  of 
using  large  values  for  the  gain  in  the  case  of  taming  by 
active  control).  In  addition  the  theory  should  be  extended 
to  include*  for  Instance*  the  possibility  of  simultaneous 
(or  quasi-simultaneous)  crossing  of  the  imaginary  axis  by 
two  pairs  of  Complex  conjugate  roots*  a  case  that  is 
excluded  here  (see  Section  2.  1).  Other  issues  thbt  should  be 
addressed  included  higher  order  nonlinear  terms*,  more 
sophisticated  aerodynamic  operators*  three  dimensional 
models*  and  adaptive  control  systems. 


APPENDIX  A 


LAGRANGXAN  EQUATIONS  OF  MOTION 
FOR  AN  AEROELASTIC  SYSTEM 


Let  the  displacement  S  of  a  point  of  an  aircraft  lie 
represented  in  terms  of  prescribed  modest  (such  as  the 
natural  modes  of  vibrations  or  finite-element  shape 
functions) 


vs  i 


(A.  1) 


inhere  t  indicates  timet 


indicates  a  set  of  material  (or 


convected)  coordinates.  The  amplitudes!  q^  ,  are  called 
generalized  Lagrangian  coordinates  of  the  system  and  satisfy 
the  Lagrangian  equations  of  motion 


A  0T  97  ^  90 


g  v  i  __  o  »  ^  -y  _  Q  (A.  2) 

M  s1i  ~  ; 

inhere  T  is  the  kinetic  energy  and  U  ■  U(q^)  is  the  elastic 
energy.  In  addition  the  generalized  forces  /gs^M|(7Minclude 
the  generalized  aerodynamic  forces  **  *  * 


Q™  s  -  A  tf.  n  .  ft* 


(A.  3> 


(inhere  A  *  '/a  V*  is  the  dynamic  pressure  and  cr  is  the 
aerodynamic  pressure  coefficients)  as  well  as  the  effect  of 
structural  dampingt  Q,  . 


Note  that 


where 


(A.  4) 


l  [h]  <*'5> 

as 

M  .  [///  f  A  *«']-[  Ill  f  MW'?] 


is  called  tha  mass  matrix.  Note  that  o I  3  |  is  time 
independent  (because  of  conservation  of  mass):  therefore* 
the  generalised  masses  are  independent  of  the  deformation. 
Therefore*  Eq.  A.  1  ensures  that  9T/3q*  *  o. 


In  addition*  assuming  that  U  is  an  analytic  function  of 
q*  and  that  q^  ■  0  corresponds  to  an  equilibrium  position 
(so  that  dU/3q.  ■  0  (i»l,...*N>  for  q.  -O*  j«l,...,N) 
yields  J 


where 


6  -[^5,1 


d- 


(A.  7) 


is  called  the  stiffness  matrix 


i* 

and 


nonlinear  terms,  which  (because  of 


O 

the 


is  of  the 


d) 


type 

.ft* 


(A.  8) 

is  the  vector  of  the 
analytic  nature  of  U) 


m-  [  i 


AJ 
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j*' 


(S) 
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+  H  Z.  £ 


For  simplicity*  it  is  assumed  that  the  aerodynamic 
forces  are  analytic  functions  of  q.  and  q  only*  so  that 

X  i 


(A.  10) 


where  g_  is  the  vector  of  the  nonlinear  aerodynamic  terms 
of  the  type  • 

2N  2N 


f-an-  { 1 1  %  - 1 


(A.  11) 
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Finally*  tha  structural-damping  forces  art  approximated 


=  B  i  + 1  '"M 


(A.  13) 


whara 


'  s  JtJ  in  *«a/  sn  % 

f  *<yj +  x^a*'"]<a' 

Combining  tha  abova  aquations  ona  obtains 
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(A.  15) 
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with 


K.1;  =  +  ^ij M 

CU\ jk.fr  )  -  +  c^l^)  +  ciflji 


(A.  17) 


(whara  b*jj  and  cf*j„  art  undarstood  to'  ba  equal  to  zero  for 
i»  or  j*  or  k  greater  than  N) 

Equation  A.  15  may  ba  recast  in  a  form  which  is 
convenient  for  tha  analysis  of  Section  II.  Setting 


^  s  *r 


Eq.  A.  15  may  be-  written  as 

la  •  .w  j.  W  a  4  A  £,  4  4  £  aT  s  ^ 
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The  above  two  aquations  may  ba  combined  to  yield 

J  k  r  A(\)  K  +£ 

where 


where  m.  and 


X  a 

3  a  | 

III] 

AfX) 

•Kij. 

f  = 

the  modes*  M.  , 
i*  then 

M  ^ 

r  •*] 

k  . [»*<**] 

are 

called  the 

(A.  21) 

coincide  with  the  natural 


(A.  22) 

i  and  thi 


natural  frequencies  respectively.  The  modes  may  be 
normalized  so  that  m,;  «  1.  In  this  case*  J  ■  I  and  Eq.  A.  20 
is  equal  to  Eq.  1.  Zn  general*  premultiplication  of  Eq.  A.  20 
by  J*'  yields  an  equation  formally  equal  to  Eq.  1. 


Zf  active  control  is  included*  the  system  is  called 
aeroservoelastic.  The  control  equations  are  usually  cast  in 
the  form  of  Eq.  A.  20  so  that  the  combined  system  may  be 
easily  recast  in  an  equation  formally  equal  to  Eq.  1. 


APPENDIX  B 


AIRFOIL  FLUTTER 

B,  %  FgT^vlilign-Bf-Profeltm 


Consider  ths  airfoil  shown  in  Figure  B.  1.  Ths  center  of 
mass*C«  is  allowed  to  move  only  in  the  vertical  direction: 
hence  the  motion  of  the  airfoil  is  limited  to  two  degrees  of 
freedom:  the  vertical  displacement!  b§,  of  the  center  of  mass 
C  (where  b  is  the  semichord  of  the  airfoil)*  and  the 
rotation  around  the  center  of  mass  C.  Let  b£  and  b*£  be  the 
Cartesian  coordinates  in  a  frame  of  reference  filed  with  the 
airfoil  (i.e. »  §  and  i»  are  convected  coordinates)  such 

that  origin  is  in  C  (center  of  mss)  and  £-axis  goes  through 
the  point  E  (point  of  attachment  of  the  springs*  i. e.  * 
elastic  axis). 

The  location  of  a  point  (  §  *  )  in  an  absolute  frame  of 

reference  (**y)  is  given  by 

X  *  +  *2,  sin*  ) 


y  *  Sinflt  4  e«i  «  +  S  ) 


(B.  1) 


The  motion  of  the 
equation  of  motion 

airfoil 

is 

governed 

by  the  Lagrange 

• 
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-  21  4 
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it 
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5Ur( 

■2L) 

HU 

du 

(B.  2) 

where  • 

indicates 

differentiation  with 

respect  to  a 

nondimentional  time 


(B.  3) 


(where  a*  is  the  speed  of  sound  of  the  undisturbed  fluid) 
whereas  T  and  U  are  the  kinetic  and  potential  (elastic) 
energy  respectively*  and  Q*  and  QA  are  such  that  the  virtual 
work  TW  is  given  as 

S  W-.  Cs  $5  +  CB.4) 

In  order  to  obtain  the  expression  for  the  kinetic 
energy*  note  that \ 

Us  x/fc  =  A,  (-  4  4i*VA  4  1  *  \ 

.  * 

nr  *  a*  y/fc  -  a*  (.  \  c*i<  -  +am& 
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Hence 
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+  [<f-£  e*i*  -  +iAtL)*  *  5]Zj  d£***L 

5i*- b*[Wjf  (i*+i1)‘*£dn.  *  ff  ?(-$£**<■ -?*««)  *£*1 

+  ikfl?dt*l]  <B6> 


or  noting  that*  by  definition  of  center  of  mass.  C  (recall 
that  £»  i^O), 


//  r  s  5  JFj  t4,irft 


(B.  7) 


one  obtains  the  desired  expression  for  the  kinetic  energy 
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where 


**  5  faa/f  f 

is  the  mass  of  the  airfoil*  whereas 


(B.  9) 


Is  b4  Jf  f  (J$‘+rlZ)d£et*L  (B.  10> 

is  the  moment  of  inertia  around  C. 

Next  consider  the  potential  energy  U.  The  motion  is 
assumed  to  be  restrained  by  two  springs  attached  at  point  E. 
For  the  first  spring  the  potential  energy  U  is  a  function  of 
the  nondimensional  vertical  displacement  £e  ■  hg/b  at  E  and 
the  second  one*  the  potential  energy  is  a  function  of  the 
rotation^.  Hence  the  potential  energy  is  given  by 


U-.  a (5,)  +  c*,(*0 


(B.  U) 


IF 


“•calling  that  tha  4-axis  goas  through  E<  so  that  «■ 
ona  obtains  ^ 


with 


Zb  =  Ut  /i  -  +  S  Slrloi  ) 


^*-S 


( B.  12) 


(B.  13) 


Naxtf  nota  that  if  p  is  tha  pressure  and  $  tha  arclanght 
along  tha  airfoil* 


£W=  c£  (~P  n  •  <T  x  )  ds 

-  _  ^  p  (n„  £  x  +  <T  if )  d  s 


Using  Eq.  1  and 


n*  *  c-cri  *  +  5/M.ai. 

11  ^  x  •  ^ ^  I1|^  fld 


(B.  14) 


(B. 13) 


(where  n$  and  na  are  tha  components  of  n  in  tha  <5*7 >  frame 
of  rafaranca)  one  obtains  (see  Eq.  B.  4) 

SVJ  ^  p  [fn-Crl*  ^  *;««)(-  Scc+rjtriol  S«) 

+  (-  ki5  *.•  -  *  4  n  t  X-  £  ***  *  «  -  *1  5/'1*  +  L 

*  -fca  ^  p[Cl\-  £  )*oL  +(-"**''**  *4*) 

s  G^ctg  +  0„  Ate  (B-  16> 

with  ds  ■  ds/b  and 

42s  a  -  fc»*  ^  jp(~nl  to** )4s  S  -  yen^  rfs 

«*  *  -fc*^  P^"S-5"l)rf*  (B.  17) 

Finally*  combining  Eqs.  B.  2*  B.  8*  B.  11  and  B.  12.  ona 
obtains 

••  a 
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J.  ®  V  _  A  (B.  18) 


with 
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V  r  I  />n  <o  (B  19) 

where  Qg  and  are  given  by  Eq,.  B.  17  and  U  is  defined  by 
Eq.  B.  11. 


B.  2  Aerodunamic  Forces 

In  order  to  complete  the  formulation*  we  need  an 
expression  for  the  pressure  distribution.  Ue  will  asssume 
that  the  air  flow  around  the  airfoil  is  supersonic  and  that 
the  pressure  is  given  by  the  piston-theory  expression  (Ref. 
If*) 

t  -(it  Jfci-  ±  )“ 

fit,  X  -Z  '  (B.  20) 

where  aM  is  the  undistured  speed  of  sound*  a'  -  cp  /c*  i»  the 
specific  heat  ratio  and  tp  is  the  normal>wash  (velocity 
component  in  the  direction  normal  to  the  surface)*  that  is 

vi/lM-  +  <**/<**  )■  ”/Oz-ng  +(*  nM 

■  00 

=  +r)  jL  (-§  vU*  + 1^*)* 

+  +  ij^^si’na+^cHe) 
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where 


(B.  22) 


For  tho  ukt  of  simplicity.  in  ths  following  wo  will 
assume  that  tho  airfoil  has  zoro— thickness  and  zero-camber. 
In  this  case.  "  ©•  Equations  B.  1  and  B.  4  remain 

essentially  the  same  (expect  for  the  fact  that  nm  0).  In 
particular.  Eq.  B.  1  reduce  to  L 

X  =  t  §  c*jeC 

y  _  ^  ^  DC  t  ^  (8*  23k 

Next,  note  that 


n^«  0 


n  «  +1  (upper  side) 

t 

■  -1  (lower  side) 

Using  these  equations.  Eq.  B.  19  reduces  to 


=  X-  C+^OC  f 

•4-y  J  a  ^ 


S  J  ~ 

-4*  J  S 

where  the  pressure  dicontinuity  is  given  by 


(B.  24) 
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the  mass  ratio  6  is  defined  as  (recall  that 

&  ^  *  .4±y~. 


(B.  26) 


(B.  27) 


and  the  integral  are  understood  from  leading  edge  to 
trailing  edge.  The  value  of  the  nondimensional  normal wash 
evaluated  on  the  upper  side  of  the  airfoil  is  given  by 


=  *  *.W+;i  C*«S  -J. 


(B.  26) 


(the  opposite  sign  holds  for  the  lower  side  of  the  airfoil). 


The  pressure  may  be  written  as 


and  (noting  that  ^  has  opposite  sign  on  the  upper  and  lower 
sides)  the  pressure  difference  is  given  by 

Ap/H.  =  *2  y  (*■$■  +  j£  M3Tfr^...)  <B.30) 

with  X  *(#+0/2- 

Combining  Eqs.  B.  28  and  B.  30  one  obtains 
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and.  substituting  in  Eq,.  B.  29. 
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^  is  the  value  of  jr  at  midchord. 
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where 


Equations  B.  32  and  B.  33  mag  be  racast  as 

Me  =  U* 


whereas 
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It  mag  ba  notad  that  in  tha  formulation  usad  bg  Smith 
and  Morino  (Raf.  19)  y»  raprasants  tha  downwash  (valocitg  in 
tha  diraction  of  tha  g-axis)  instaad  of  tha  normalwash 
(valocitg  in  tha  diraction  of  tha  v^-axis);  hanca  tha 
axprassion  for  f>  consistant  with  tha  formulation  of  Ref.  19 
is 


_  J_  5  U  CoSoi 


(B.  38) 


(saa  Eq.  C6  of  Raf.  19) 

Using  this  axprassion  one  obtains  that  e^  and  <  are 
still  given  bg  Eq.  B.  39  with  # 
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raspactivalg  (in  agreement  with  Eqs  C12c  and  C12d  of  Ref. 
19) 


B,  3.  ElajLils  Efluai 

For  tha  sake  of  simplicitg  tha  problem  is  limited  to  tha 
case  in  which  tha  potential  energies  U^  and  U^  are  even 
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functions  of  and  oC  respectively  (this  implies  for 

instance  that  the  axial  spring  acts  symmetrically  in 
compression  and  tension)#  and  are  given  in  the  forms 

U  =  '/a  U*  4«/f  Kj  ^  +... 

U*  ^  Vi  Kk  +y4  K*  **  +*•  • 

Using  Eqs.  B.  12  and  B.  19  one  obtains 
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These  equations  may  be  rewritten  as 
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where 
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in  agreement  with  Ref.  19  (Eqs.  C12a  and  C12b). 
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Combining  Eqs.  B.  18,  B.  39,  and  B.  47  on*  obtains 

M/  +5/  +  (K  +A  E)  y  s  2 
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P..3.  fiotfif guattm 

In  order  to  achieve  flutter-taming,  the  system  may  be 
modified  for  instance  by  changing  the  design  to  obtain 
different  values  of  the  nonlinear  stiffness  coefficients 
kr  and  k lc  .  If  the  structure  is  already  built  (such  as  in 
the  case  of  existing  airplane)  one  could  add  a  nonlinear 
damper,  causing  for  instance  a  moment  proportional  to  it3 , 


i.  *.  , 

modif ing 

the  expression  for  gA  (Eq.  B.  97)  as 

(B.  98) 

with 

(B.  99) 

A 

more 

ambitious  approach  (an  active 

control 

modification)  may  be  obtained  by  "Using  a  control  surface. 
Indicating  with  the  deflection  of  the  control  surface,  the 
normalwash  is  given  by  (see  Eq.  B.  28) 
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where  J„is  tha  location  of  the  hinge  line.  For  the  analysis 
presented  here  (in  which  9  is  of  the  order  of  %* or  ot*  )  the 
above  expression  nay  be  approximated  by 

*i/<  «*  +  (5  -5*  -'/» 

-WCS-Sh)  e  -*<5-5»)  e  +•••  <B  «> 

where  H(u)  is  the  Heaviside  function  and  R(u)  >uH(u)  is  the 
ramp  function.  The  last  two  terms  yield  an  additional 
pressure  distribution*  p  (due  to  the  control  surface 
deflection  &  )#  equal  to 


f»*  -2?  [M  Hft-f,,)®  +-R(5-5h)6J 

This  yields  an  Increase  in  qf  and  q^  equal  to 
A  fa  s  €  (rtSf'e  +  $“>  e  ) 
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where 
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and  >▼«  y* 

where  d  ■  (  %n-  §h  )/2  is  the  ratio  between  the  chord  of  the 
control  surface  and  the  chord  of  the  airfoil. 

The  dynamics  of  the  control  surface  is  assumed  to  be  to 
irreversible  (i.e.  #  unaffected  by  the  aerodynamic  forces) 
and  governed  by 


J9  +  re  +  ke  = 


(B.  67) 


where. the  feedback  term  is  proportional  for  instance 

to  the  vertical  displacement  of  the  center  of  mass: 
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The  resulting  equations  ere 
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with  g  given  by  Eq.  B.  62  end  g^  given  by  Eq.  B.  70 
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It  should  be  emphasized  that  the  three  modifications 
introduced  in  this  section  (a.  modify  k's  and/or  k^  »b.  add 
a  damper  producing  a  moment  proportional  to  &s  < and  c. 
nonlinear  feedback  used  in  the  dynamic  equation  for 
deflection  of  the  control  surface)  are  only  examples.  Many 
other  modif icationst  equally  vail  d«  can  be  devised  (such  as 
a  damper  creating  a  vertical  ~force  proportional  to  ). 
However  these  examples  should  be  sufficient  to  give  a  good 
idea  of  the  flexibility  available  in  implementing  the  idea 
of  flutter  taming. 
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